Bulk-edge coupling in the non-abelian = 5/2 quantum Hall interferometer 
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Recent schemes for experimentally probing non-abelian statistics in the quantum Hall effect are based on 
geometries where current-carrying quasiparticles flow along edges that encircle bulk quasiparticles, which are 
localized. Here we consider one such scheme, the Fabry-Perot interferometer, and analyze how its interference 
patterns are affected by a coupling that allows tunneling of neutral Majorana fermions between the bulk and 
edge. While at weak coupling this tunneling degrades the interference signal, we find that at strong coupling, 
the bulk quasiparticle becomes essentially absorbed by the edge and the intereference signal is fully restored. 



PACS numbers: 73.43.Cd, 73.43.Jn 

Recently, interference experiments were proposed as a way 
to examine the non-abelian nature of quasiparticles in the 
v = 5/2 quantum Hall state iH |2, d |3l. The most dramatic 
signature of non-abelian statistics is expected to be seen in 
the interference of back-scattering amplitudes from two con- 
strictions in a long Hall bar. (See inset in Fig. 1.) The two 
constrictions enclose a "cell", whose area may be varied by 
means of a side-gate. The bulk is assumed to host a number 
Nqp of localized quasiparticles, that do not take part in elec- 
tronic transport, and have no tunnel coupling to the edge. In 
the limit of weak back-scattering, when Nqp is even the two 
back-scattering amplitudes interfere coherently, while when 
Nqp is odd they are incoherent, and thus do not interfere. In 
the former case, the back-scattered current oscillates with the 
area of the cell, while in the latter case it does not. This dif- 
ference reflects the non-abelian nature of the quasiparticles. 

The theoretical analysis makes a sharp distinction between 
bulk and edge. In a real system, however, some degree of 
coupling between the edge and quasiparticles localized in the 
bulk is unavoidable. Since both the 5/2 edge and the quasipar- 
ticles consist of both neutral Majorana fermionic and charged 
bosonic degrees of freedom, several types of edge to bulk cou- 
pling are possible. We expect that at low energies, tunneling 
that involves a charge will generally be suppressed due to the 
Coulomb energy. Thus in this work we will focus on tunnel- 
ing of the neutral Majorana mode from the bulk to the edge, 
and on the resulting effect on the interference. 

The system we consider IH H S 0] is a Hall bar lying par- 
allel to the X-axis (See Fig. 1). Two constrictions are located 
at X = —b and x ~ b. We focus on a simple case where there 
are two quasiparticles, Nqp = 2, localized at a; = 0, between 
the two constrictions, with one of the quasiparticles coupled to 
the upper edge and the other coupled to the lower edge. Ref. 
101 considers the case of Nqp = 1 in the weak tunneling limit. 

When the two localized quasiparticles are decoupled from 
the edge they form a two level system, and the ground state 
is doubly degenerate. The interference patterns that are seen 
in the two respective ground states are mutually shifted by 
a phase tt. Then, at temperature T = the magnitude of 
the interference term depends on the ratio of two time scales. 



one determined by the voltage ty = h/e*V, and the other 
being the time associated with motion between the two con- 
strictions tf, = 2b /v where v is a characteristic edge mode 
velocity. When analyzing the effect of bulk-edge coupling we 
will focus on the case of low voltage, ty ^ h, where the 
interference is most clearly seen. 

We start with a qualitative description of our results. In 
the absence of edge-bulk coupling the system cannot switch 
from one ground state to another. Thus if it is prepared 
in one ground state, repetitive measurements of the interfer- 
ence would show the same interference pattern. However, 
when the interference term is averaged over the two possi- 
ble ground states, e.g., by measuring the interference with 
a random choice of the initial ground state, the average is 
zero. When the coupling of the bulk two-level system to 
the edge is turned on, the average value of the interference 
term becomes non-zero, and the correlation function between 
consecutive measurements is strongly modified. Denoting 
the coupling strengths between the localized Majorana par- 
ticles and their respective edges by A„ and A^, [defined pre- 
cisely in Eq. ([T]i below], we obtain corresponding time scales 
t\u{d) = ('^^u(d)/2«™)~^ where u,„ is the velocity of the 
Majorana modes on the edges. In the limit of weak coupling, 
where txu{d) ^ ty, we may use a perturbation analysis, and 
we find that the average value of the interference is propor- 
tional to {tvy^'^itv ^og^ltx/tv], where we have assumed 
that t\u and txd are comparable in magnitude, and t\ is their 
geometric mean. As the coupling is increased, or as the volt- 
age is lowered, the perturbative analysis breaks down. We 
then carry out a numerical analysis, which suggests that in the 
limit tv /t\ — > oo the full magnitude of the interference term 
is retrieved. In effect, the two bulk quasiparticles become then 
a part of the edge, and Nqp reduces from two to zero. (We find 
a similar effect for a single quasiparticle strongly coupled to 
an edge.) In contrast to the build-up of the average interfer- 
ence term as the coupling gets stronger, the fluctuating part of 
the interference pattern is weakened by the coupling, and its 
characteristic correlation time becomes tx, which decreases 
with increased coupling. 

For the derivation of these results we will follow several 
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steps. After the introduction of the relevant Lagrangians, 
we derive the operator that describes quasiparticle tunneling 
across the two constrictions in the presence of the two local- 
ized bulk quasiparticles. We find it useful to represent this 
operator in two forms, a local form using the a operator of 
the Ising Conformal Field Theory (CFT) that describes the 
V = 5/2 edge, and a non-local form in terms of the Ma- 
jorana fermions that propagate along the edge. Within the 
non-local form, we show that the tunneling operator is pro- 
portional to a "parity operator" that measures the parity of the 
number of electrons encircled by a back-scattered quasiparti- 
cle as it moves from x — —oo along one edge, through the 
back-scattering at the constriction, back to a; = — oo along 
the other edge. Next, we perturbatively analyze the weak cou- 
pling limit, and finally we numerically analyze the strong cou- 
pling limit. 

In the absence of any coupling to bulk quasiparticles the up- 
per (w) and lower (d) edges of the = 5/2 state are described 
by two charged boson fields 4)u{x), 4>d{x) and a neutral Majo- 
rana fermion field. The Lagrangian for the boson field on each 
edge is that of a chiral Luttinger liquid, characterized by a ve- 
locity -iiVc- The Lagrangian density for the Majorana fermion 
field is £^ = -^4>^{x)[dt — v^^dx]i^'' {x) with r taking the 
values u and d for the upper and lower edges. For simplicity 
we set the velocities of the Majorana edge modes to be equal 
and opposite v!J^ — ^"^m ~ '^m- Furthermore, we set Vm — ^ 
when no confusion results. The Majorana Lagrangian can also 
be thought of as the Lagrangian of an Ising CFT 1911 . 

Each of the two localized bulk quasiparticles carries a zero 
mode, described by a localized Majorana operator. We de- 
note the two bulk Majorana operators by 7^ , 7^, with the sub- 
script indicating the edge to which the quasiparticle couples. 
The two-dimensional Hilbert space created by the two Majo- 
rana modes is spanned by the two eigenvectors of the operator 

To examine the effects of bulk-edge coupUng we couple 
to the upper edge and 7^ to the lower edge, both at a; = 0. The 
Lagrangian density for this coupling is 

£b_e = I [\uV{x)lu + \d^''{x)ld] S{x) . (1) 

The Lagrangian £"j + + £f,_e introduces the time scales 
txu{d) defined above. The bulk-edge coupling mixes the states 
with eigenvalues ±1 of i^uld- Roughly speaking, t\ is the 
time in which a state with a particular value of i^uld decays 
to a mixture of the two eigenvalues. 

The operator that tunnels a quasiparticle across a constric- 
tion may be expressed in a local form through the cr,i , cjd op- 
erators of the Ising CFT that describes the v = b/2 upper and 
lower edges iQ]. The tunnelling operator is Htun ^ T + T+, 
where 

T = e*^* ^* [t^lClMl + mCRMBhuld] (2) 

transfers a quasiparticle from the lower to the upper edge 
through the left (L) and right (i?) constrictions respectively, 
and its hermitian conjugate similarly transfers a quasi- 
particle from the upper to the lower edge. Here, V is the 



voltage difference between the two edges, e* = e/4 is the 
quaisparticle charge. Correspondingly, the current operator 
is given by J = ^(T — r+). The operators C/^j^) = 
gi(0„(T&)-0d(Tb))/V8 tjje charge part of the tunneling op- 
erator, operating on the charge mode. The Aharonov-Bohm 
phase is absorbed into the relative phase between the tunnel- 
ing coefficients r/^ j^. The neutral parts of the tunneling oper- 
ators areWi = ct„(-6)ct<j(-&) and TV}? = cr„(fo)CT<j(6). For 
the present purpose, the a operators are defined through their 
operation on the Majorana fermion fields as 1^ 

(Tr{xo)il^r{y)crr{xo) = -sgn(a;o - y)i'r{y) (3) 

with r ~ u,d. The factor of ju'^d in the second term of Eq. (|2| 
is included to account for the wrapping of a tunneling quasi- 
particle at position x ^ b around the two localized quasiparti- 
cles. This factor is responsible for the tt phase shift between 
the interference patterns corresponding to the two eigenvec- 
tors of 7„7rf. 

The neutral mode part of the tunneling operators may 
also be expressed in a non-local form through the Majorana 
fermions along the two edges in a way which we find to be 
both illuminating and useful. This approach is based on the 
description of the u = 5/2 state as a p-wave superconductor 
of composite fermions [1]. Within this description the bulk is 
a superconductor, with the localized quasiparticles being vor- 
tices in that superconductor. A tunneling of a quasiparticle 
from one edge to another at position xq involves a tunneling 
of a vortex, and that introduces a twist into the phase of the 
order parameter: for all points in the region x < xq, the phase 
is shifted by 2tt, while for all points in the region x > xq the 
phase is unaffected by the vortex motion (up to an unimpor- 
tant global gauge redefinition). To implement this shift of the 
phase, we recognize that the phase field is canonically conju- 
gate to the Cooper-pair density field, which at zero tempera- 
ture is just half the electron density field. The operator that 
implements the required shift in the phase is then 

7-1/ \ r ^ dr p(r) 

P(— oo,xo)=e ■'^^==0 . (4) 

Since the operator dr p(r) has only integer eigenval- 

ues, the operator P(— oo,a;o) is nothing but a Parity Oper- 
ator which measures the parity of the number of electrons 
to the left ifioll of xq. Eq. (|2]i can thus be rewritten as T = 

e'^*^* [77lCl-P(-oo, -b) + 77flCRP(-oo, b)] as we shall see 
below. 

Since the bulk of the system is gapped, and since all parti- 
cles in the superconducting ground state are paired, the parity 
operator only has contributions from localized neutral modes 
and from the neutral mode along the edge. The operator 
ijuld in the second term of Eq. (|2]i precisely counts the par- 
ity of the number of fermions in the localized bulk quasi- 
particles to the left of X ~ b. Counting the fermions along 
the edge is a bit more complicated but is achieved by con- 
structing a complex Fermi field ipe{x) = ipuix) — iipdix) and 
iplix) — 4'u{x) + iipd{x), such that the edge contribution to 
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the parity operator is 



(5) 



It is easy to see that Eq. ^ holds when the operators (Jr{xo) 
are replaced by Pedge- The eigenvalues of the latter are 
±1, since the eigenvalues of dxipl{x)'ipe{x) are integers. 
The application of either ipdiu) or ipu{y) on an eigenstate 
of dx^l{x)ipeix) changes the eigenvalue by ±6{xq — 
y), and hence ([3]). Altogether, then, we have AfLiR) = 

Pedge{-00,Tb). 

To calculate the current-voltage characteristics in the 
weak back-scattering limit, we use standard pertur- 
bation theory in the tunneling strength to yield / = 
— i j'^^dt {[J{0), Htunit)])- With some algebra, the interfer- 
ence term that results is 



Re 



dt( 



'vt 



X ([C+WAAi(0,Cfl(0)AAflX0)7«(0)7d(0)]) .(6) 

For —e*V > 0, only the first term of the commuta- 
tor, with t-dependent operators to the left, will contribute 
to the integral. The correlator of the charged operators 
(the C's) and that of the neutral operators (the JV's and 
7's) factorize. The correlator of the charged operator is 
{C+{t)CR{0)) = [S + t{v,t ~ 26)]-i/8[,5 + i{vct + 26)]-i/8, 
where (5 is a short-distance cutoff. The bulk-edge cou- 
pling affects only the neutral correlator, to be denoted by 
I^/{t) = (AAL(t)7Vfl(0)7„(0)7d(0)) which is just the parity- 
parity correlator {P{—oo,—b;t)P{—oo,b;t = 0)). In the 
absence of edge-bulk coupling, this correlator breaks into 
a product (7VL(t)7Vj?,(0))o (7«7rf)o- We denote XaA^"^ = 
{AfL{t)AfBXO))o which has the value [S + i{t ~ 2b)]-^/^[S + 
i{t + 26)]-!/*^ 11 H. The correlator (7„7d)o is ±i, de- 
pending which ground state is considered. For either ground 
state, integration of these two expressions in Eq. ^ leads to 
an interference term of the same visibility as in the absence 
of any bulk quasiparticles ll ill . Since we are interested in 
the effect of the bulk-edge coupling on the visibility of the 
interference, we find it useful to define a reduction factor 
R{t)=lM{t)/lM'^'\t). 

We now turn to analyze the reduction factor in various 
regimes of bulk-edge coupling. Generally, the two edge theo- 
ries (u, d) factorize and we can write the correlator Tj^{t) = 
Tu{t)Td{t). In the limit of weak coupling, we may use per- 
turbation theory. We expand the time evolution operator to 
lowest order in A. The perturbed correlators can be written as 
correlators in an unperturbed theory 

luit) = (7) 
Xufdt' {Tc7.,i~b,t)a^{b,0)^u{0,t'))o (T7„(0)7„(t'))o 



where the time integration contour starts at — oo goes up to 
t across the real axis then back to — oo and T represents the 
appropriate (Keldysh) time ordering of operators. 



The correlators of the type {aaip)o are well known 
from conformal field theory Jg]]: {a-{zi)a{z2)ip{z3))o ~ 
^(212)^/^ (^23213)""^''^ where z,y = z.^ - Zj and z = x + ir 
in imaginary time, which then needs to be continued back to 
real time. Substituting this correlator in Eq. (|7]i and noting 
that at the unperturbed level (7M(0)7u(t))o = 1 we find the 
time integral to be logarithmically divergent. However, when 
the correlator {ju{0)'yu{t)) is itself calculated in perturbation 
theory, it is found to decay at a time scale of order t\. Thus 
this correlator provides a natural cutoff for the time integra- 
tion. Evaluating the integrals with the cutoff yields (in the 
limit of small tf,) that the leading contribution of the upper 
edge to the parity correlator (which is independent of the de- 
tails of the cutoff) is 

luit) = (5 + iif/«{A„%/2Hlog(|t|AAj-7rsgnt]}. (8) 

When we consider coupling of impurities to both edges, 
we obtain a similar expression for Td- The reduction factor 
defined above is then 



i?(t) = 2A„A,tlog(|t|AAjlog(|t|AAj 



(9) 



Including the contributions from the charge modes and 
in Eq. results in an interference current proportional 
to XuXdV^"^^^ log{tx^)\og{tx^). Interestingly, in the case 
where there is only a single bulk quasiparticle, coupled to 
just one edge, the corresponding logarithmic factor disappears 
from the interference current, which, in the weak coupling 
limit, is proportional to XV^^, as was shown in Ref. ll7i ll2ll . 



0.96 



0.92 



0.88 



0.84 




FIG. 1: Reduction R{t) of the imaginary time parity-parity correla- 
tion function due to coupling between edge and localized modes. In- 
set shows the interferometer with two localized Majorana modes -ju, 
"fd coupled to upper and lower edge, respectively. Data were numeri- 
cally obtained for an interferometer of size b = tx/'l and total system 
size L + b — 10.5tA, where t^^ is the characteristic decay rate of 
localized modes. Data points are shown for a/tx ~ 0.005 (full line 
with open circles), a/tx = 0.0025 (plus signs), and a/tx ~ 0.00125 
(asterisks), where a denotes the lattice constant. 

In order to analyze the strong coupling limit with either ti, 
or ty of the order of tx, we numerically study a lattice ver- 
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sion of our model lll3ll . We start with a tight-binding Hamil- 
tonian for one-dimensional complex fermions without the lo- 



a is the lattice constant, and the operators obey the usual anti- 



calized modes: H = —VmO, ^^Ac\^-iCj + h.c). Here, 



commutation relations {cj, Cj} = 5ij. We study the model at 
half filling with a Fermi wave vector kp = 7r/2. The fermions 
created by can be decomposed into two Majorana species 



7j 



Cj + e c] 



We define continuum fermions, using 7^ alone, via 

= -^5](±1F7(j-/)7/ , (10) 

where /(j) is a Gaussian with a width large compared to the 
lattice spacing and subject to the normalization f{i) = 
-^/tt. Using this mapping one can now include the coupling to 
the localized modes as in Eq. ([T]). 

The parity operator for a set of lattice sites {j} can be writ- 
ten as the product over sites of operators 2cjcj — 1 = ijj^j- 
The parity operator for the localized modes has a similar form. 
The expectation value of any such product, at the same or dif- 
ferent times, can be evaluated, using Wick's theorem, as the 
Pfaffian of a matrix whose elements are the pair correlation 
functions of operators on different sites, including the local- 
ized modes where appropriate. As the Hamiltonian is a sum 
H = + H^, the 7 species contributes a factor to the par- 
ity expectation value which is the same whether the localized 
modes are present or not. Thus we may ignore the 7 modes 
in evaluating the reduction factor In this way, the lattice form 
for the edge parity for a region [xi , X2] becomes 



Pedgeixi,X2) 



n 



7j 



(11) 



Xi<Xj<X2 



To calculate the reduction factor at time t we will 
need to evaluate expectation values of products like 

{Pedge{-L,-h]t)Pedge{-L,b]0)i^uld) , whcrC -L <C -|t| 

is a point far to the left of the origin. In the absence of bulk- 
edge coupling, the pair correlation function for two lattice 
points at equal times is 



(«7j7fe)o = -- — 



TT J 



(12) 



The two terms in brackets arise from the right-moving up- 
per edge and left-moving lower edge, respectively. With non- 
zero bulk-edge coupling, the correlation functions (i'jj'jk) 
and {ijjju{d}) can be calculated analytically in the contin- 
uum limit (i.e., for points not too close to the origin) for the 
same and for different times. A full lattice calculation can be 
carried out numerically, but it is time-consuming for large lat- 
tices. We have found that the short-distance errors introduced 
by using continuum correlation functions only lead to an error 
in the Pfaffian by a factor c\ that depends ona/t\ but is inde- 
pendent of the interferometer size b and the time difference t. 



The numerical results presented in Fig. 1 were obtained using 
the continuum correlation functions and corrected by the fac- 
tor cx . The numerical value of cx can be obtained most easily 
by considering equal time correlation functions. As a check, 
we note that results for different values of a/t\ obtained with 
this numerical technique deviate less than 0.3 % from each 
other, as shown in Fig. 1. 

Fig.[T]displays the imaginary time reduction factor R{t) for 
intermediate and strong coupling, for an interferometer size 
h — <a/4. R{t) is related to the real time reduction factor 
via the analytic contiuation R{t) = R{t it + 6). The re- 
duction factor R{t) monotonically increases with increasing 
time. At T « t\, there is a crossover from parity reduction 
determined by the interferometer size 6 = f a /4 to parity re- 
duction determined by the time. At large times R{t) seems to 
saturate at a value of one, implying that its analytic continua- 
tion R{t) saturates near one as well. Note that when i?(t) = 1 
the visibility of the interference is the same as it would have 
been in the absence of the two bulk quasiparticles. Similar 
results are expected if we have one strongly coupled localized 
mode inside the interferometer path, and a second localized 
mode, of arbitrary coupling, outside the interferometer We 
attribute the re-emergence of the interference as the bulk-edge 
coupling gets strong to the correlations that develop between 
the occupation of the fermionic mode associated with the two 
quasi-particles and the occupation of the region of the edge 
at a distance Vmt\ from the coupling point. Each of these 
occupations strongly fluctuates due to the coupling, but their 
fluctuations are strongly correlated. 

In conclusion, we have found that when the coupling be- 
tween the Majorana mode associated with a localized e* = 
e/4 charged quasiparticle and an adjacent edge is sufficiently 
strong, so that the characteristic tunneling time t\ is short 
compared to the time scale ty = he* /V set by the voltage, 
it appears as if the localized quasiparticle has become part 
of the edge. Specifically, for an interference path enclosing 
the quasiparticle, the interference visibility should have essen- 
tially the same strength as if the quasiparticle were not there. 
For weak coupling, the time-averaged interference intensity 
is reduced, by a factor which is oc {tv/tx)\og'^{tx/tv) in 
the case where there are two localized quasiparticles inside 
the loop, coupled respectively to the two edges with similar 
strength. 
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